Ellipsoidal Coulomb Crystals in a Linear Radiofrequency Trap 



U. Frohlich, B. Roth, and S. Schiller 

Institut fiir Experimentalphysik, Heinrich-Heine-Universitdt Diisseldorf, 40225 Diisseldorf, Germany 



A static quadrupole potential breaks the cylindrical symmetry of the eftective potential of a 
linear rf trap. For a one-component fluid plasma at low temperature, the resulting equilibrium 
charge distribution is predicted to be an ellipsoid. We have produced laser-cooled Be^ ellipsoidal 
ion crystals and found good agreement between their shapes and the cold fluid prediction. In two- 
species mixtures, containing Be"*" and sympathetically cooled ions of lower mass, a sufficiently strong 
static quadrupole potential produces a spatial separation of the species. 



One-component plasmas have attracted significant at- 
tention in the past, since they represent simple multi- 
particle systems that can be studied under a variety 
of conditions with a high degree of experimental con- 
trol. Detailed theoretical analysis, both analytical and by 
molecular dynamics (MD) simulations, is possible 0, Q]- 
To overcome the Coulomb repulsion, the plasmas are con- 
fined in Penning- or Paul- type traps. Their temperature 
can be varied over many orders of magnitude. In par- 
ticular, they can be efficiently cooled to the mK range 
by laser cooling. Strong cooling results in phase transi- 
tions to a crystalline state, whose occurrence is described 
by the interaction parameter F = Q'^ /AireoakBT, the ra- 
tio between average nearest-neighbor Coulomb energy (a: 
average particle spacing, Q: particle charge) and thermal 
energy. MD simulations on infinite systems have shown 
that the plasma becomes fluid, i.e. exhibits spatial corre- 
lations for F > 2, without going through a discontinuous 
gas-fluid phase transition |21 ■ For F > 170 0] a phase 
transition to a crystal occurs. Coulomb crystals in hy- 
perbolic Q and linear Paul (rf ) traps 0, |^ have become 
of great importance in quantum optics, where they can be 
used to implement quantum gates or quantum memories 
, and can serve as systems forprecision measurements 
on atomic or molecular ions 0, |^ . 

The spatial distribution of a trapped one-component 
gaseous plasma differs signiflcantly from the fluid and 
crystalline state. This is described by the Debye length 
— {ksTeo/nQ'^Y/'^ which is the distance at which 
interactions between individual particles overcome col- 
lective effects 0- Here, n is the particle density of the 
plasma. When A^i is much larger than the spatial extent 
of the plasma, the Coulomb interaction is negligible, and 
the density is a local function of the trap potential only. 
For a harmonic potential, the density has a Gaussian de- 
pendence on the coordinates, with extensions inversely 
proportional to the respective trap potential curvatures. 
When Ad is comparable to or smaller than the spatial 
extent of the plasma (as is the case for the plasmas pre- 
sented in this paper), space charge becomes important. 
The density is then a local function of both trap po- 
tential and space charge potential. Since the latter de- 
pends on the density distribution over the whole space. 



a self-consistent density distribution arises, whose shape 
exhibits a non-trivial dependence on the trap potential 
curvatures. For cylindrical symmetry, such plasmas are 
spheroids (ellipses of revolution) and have been studied 
both in linear rf and Penning traps [lol . In absence 
of cylindrical symmetry the shape has been predicted to 
be that of an ellipsoid Ellipsoidal plasmas have al- 
ready been observed and studied in Penning traps (see 
[l^ and references therein) . In this work we describe the 
observation and characterization of cold ellipsoidal plas- 
mas in a linear rf trap. In particular, we have obtained 
ellipsoidal crystals containing two different ion species. 
The linear rf trap used in this experiment consists of 
four rods of radius r' (inset in Figl^J, each divided into 
three electrically isolated segments. A radio frequency 
voltage Vrf sin Q.t and a static voltage Vdc are ap- 
plied to the rods in a quadrupolar configuration. Con- 
finement along the trap axis (z-axis) is achieved by 
raising the two end segments of each rod by a static 
voltage Vec- When the Mathieu stability parameter 
q = 2QVrf I Tnar^^"^ <C 1, the independent motion of 
the trapped ions (mass m) can be adequately described 
by the motion in an harmonic effective potential (pseu- 
dopotential) C/trap(r) = 'm^ujfxf /2 [ij, and an addi- 
tional jitter motion at the radio frequency f2, the so- 
called micromotion, tq is the minimum distance from 
the electrode surfaces to the trap axis. The axial (uJz) 
and transverse {cOx and tUy) frequencies of the effective 
trap potential are given by — 'inQVEC l^n, '^t y — 
Q'^Vlpl'^^rn^^^rl) - ujI/2 ± QVoc/mrl, where k is a 
constant determined by the trap geometry. For vanish- 
ing static voltage Vdc, and LOy are degenerate. With 
the application of a static voltage Vdc Oj a static 
quadrupole potential is added to the effective trap poten- 
tial and the cylindrical symmetry is broken. The trans- 
verse trap frequencies ujx and Wy increase and decrease, 
respectively, until Wy vanishes for a sufficiently large ap- 
plied voltage Vdc, which implies that the ion motion 
along the y-direction becomes unstable. 
In the fluid phase, a trapped plasma in thermal equi- 
librium at a given temperature T may be regarded as 
a macroscopic charged fluid with number density nf(v) 
This description leads to an expression for e.g. 
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the equilibrium shape of the trapped plasma. A particu- 
larly simple analytical treatment is possible in the limit 
of an ultracold plasma, when T 0. In this case, the 
equilibrium number density nf(r) is determined by the 
condition of no net force on any plasma region. How- 
ever, at very low temperatures, the plasma does not re- 
main fluid, but crystallizes with the appearance of shells 
and a fairly complicate order. While the fluid description 
may be expected to become inaccurate in this case, it has 
been shown to remain applicable as far as the shape of the 
outer boundary is concerned and as long as the crystals 
are sufficiently large. For linear rf traps, this comparison 
has so far only been performed for the case of cylindrical 
symmetry {Vdc = 0) [lO]. 

In the effective potential approximation, the absence of 
net force implies a balance between the electric field due 
to the space charge potential 4>f{r) and the trap force 
due to the effective potential: Q(j)f{r) + Utrapi'r) — const. 
The number density follows directly from Poisson's equa- 
tion, n/(r) = (eo/(3^)A[/trap(r). For the harmonic ef- 
fective potential of the linear rf trap with the above 
trap frequencies the number density is explicitly given by 
no = eoV^p/mVP'r'^, which is constant within the fluid 
and independent of the static voltages Vdc and Vec ■ 
The outer shape of the zero temperature charged fluid 
with constant number density is an ellipsoid with princi- 
pal axes RxjRy and L in the x-, y- and z-direction, re- 
spectively [131 • With the additional boundary condition 
that the electric space charge potential ^/(r) vanishes at 
infinity, the potential inside the fluid is given by 

<Pf{v) = ^[A.{Rl - x^) + Ay{Rl - y')+ 

where the dimensionless functions Ax , A, and Az depend 
on the principal axes Rx, Ry and L Force balance 
implies that the Ai satisfy ujf — Q^noAi/2meo. The 
ratios of the principal axes Rx/L and Ry/L can thus be 
calculated by solving the set of equations 

/uxV _ Ax{Rx/L,Ry/L) 
UJ Az{Rx/L,Ry/L) 

fu^y __ Ay{Rx/L,Ry/L) 

\ujj Az{Rx/L,Ry/L) ■ 

For Vdc = 0, and Uy are equal, and the equilibrium 
shape of the zero temperature charged fluid is a spheroid 
with radius R := R^ = Ry and half length L. The two 
equations l(2Jl then reduce to a single equation for the 
aspect ratio R/L. Hornekaer et al. [lO|, [llj have found 
good agreement between this theory and experimental 
aspect ratios of a large variety of spheroidal '*°Ca'^ ion 
crystals, obtained by varying the end-cap voltage Vec 
and radio frequency amplitude Vrf- 
We have tested the predictions of the charged fluid 



model in a fully anisotropic effective trap potential using 
a laser-cooled ^Be^ ensemble. The linear trap properties 
were tq = 4.32 mm, a radio frequency n/2TT = 14.2 MHz 
with an amplitude V^e = 380 V. This resulted in a 
small Mathieu stability parameter of q = 0.055, implying 
that the micromotion was relatively small and that the 
effective potential description is appropriate. For axial 
confinement, Vec = 4.5 V was applied to the trap end 
segments, giving rise to k = 3.0 • lO^'^/mm^, an axial 
frequency of lUz — 2tt ■ 85 kHz, and a transverse frequency 
of LUr = 2tt ■ 268 kHz. Experimentally, both values 
were obtained from a measurement of the transverse 
frequency as a function of Vec, in absence of Vdc- This 
was done by external excitation of the radial motion 
of gas phase ^Be^ ions in the harmonic trap potential, 
which was detected by a drop of the fluorescence signal 
from the laser cooled ^Be"*" ions [l7j . 
When the effective trap potential is made anisotropic 
by application of Vdc the predicted instability limit 
is 4.96 V. We observe partial particle loss when Vdc 
exceeds 4.2 V and total loss at Vdc = 4.9 V. To prevent 
particle loss, we limited Vdc to a maximum of 4.2 V, 
corresponding to an increase of the transverse trap 
frequency lUx/'^tt from 268 kHz to 365 kHz and a decrease 
of ujy/2Tr from 268 kHz to 104 kHz. Experimentally, the 
two transverse frequencies Wx and ujy can be measured 
by secular excitation. 

The trap was loaded with ^Be"*" ions by evaporating 
beryllium atoms from an oven and ionizing them in 
the trap center by electron impact. The trapped ^Be"*" 
ions, initially forming a hot plasma cloud of « 1000 K, 
were laser cooled by laser radiation at 313 nm until they 
finally underwent a phase transition to a crystalline 
state with a temperature of a few mK. A description 
of the all-solid-state laser system is given in To 
image the ^Be"*" ion crystals, a CCD camera was placed 
transverse to the trap axis. FigQ shows an ion crystal 
containing ~ 2.0 x 10'^ ^Be^, at different values of Vdc- 
The estimate of the ion number in Fig^ (a) is obtained 
from molecular dynamics (MD) simulations in which 
the observed structure (especially the number of shells) 
is reproduced . For the test of the calibration of the 
CCD optics magnification as well as the determination 
of absolute dimensions of the ion plasmas, we also use 
the MD simulations. As a check, multiplying the volume 
47ri?45L/3 of the crystal shown in Fign](a) and the cold 
fluid model density ng, we obtain the value ~ 2.0 x 10"^. 
This agrees well with the MD results and implies that 
the model is applicable for large crystals. 
In the outer region of the crystal shown in Fig^the ions 
are arranged in concentric shells, while the core appears 
to be fluid. The size of the liquid core appears to increase 
with applied static quadrupole potential. A possible 
explanation for this behavior is the increased thermal 
motion of the ions in the core, which was reproduced 
by MD simulations and is a subject of ongoing studies. 




FIG. 1: CCD images of a large ion crystal containing ~ 2.0 x 
10^ ^Be+, taken perpendicular to the 2-axis and at 45° to the 
X- and j/-axes, for different values of Vdc- Dashed lines: fits of 
ellipses with principal axes 7?45 and L to the outermost shells, 
(a): Vdc = 0, the effective trap potential has cylindrical 
symmetry. The ion crystal is a prolate spheroid with radius 
Rx ~ Ry ~ Rib and half length L. The aspect ratio R45/L = 
0.178. (b)-(e): Vdc 7^ 0, the cylindrical symmetry is broken 
and the crystal is an ellipsoid with principal axes Rx < i?45, 
Ry > -R45 and L. With increasing Vdc the crystal expands 
in the j/-direction and compresses in the x- and z-directions, 
while maintaining constant volume. Vdc is set to 1.8 (b), 
2.8 (c), 3.6 (d) and 4.2 V (e), leading to an aspect ratio of 
R45/L = 0.193 (b), 0.235 (c), 0.309 (d), and 0.413 (e). For a 
definition of Rx, Ry and R4S see inset of Fig|5| 



In contrast to the observation of partially crystallized 
large plasmas, we found that under similar laser cooling 
parameters small plasmas containing only a few hundred 
^Be"'' ions always crystallized completely. The left 
central part of the ion crystal contains a dark region, 
which consists of sympathetically cooled ions originating 
from the residual gas and having a mass smaller than 
that of ^Be^, pushed to one side by radiation pressure 

The ellipsoidal deformation is a reversible process if the 
maximum value of Vdc is kept within the range stated 
above: after turning Vdc off, the principal axes i?45 and 
L of the crystal returned to the initial values, indicating 
that no ions were lost during deformation. The outer 
boundaries of the crystals in FigQlcan be well described 
by ellipses with principal axes i?45 and L. Since the 
CCD camera takes a projection along an axis at 45° with 
respect to the x- and y-axes, the principal axis i?45 of 
each ellipse is related to the principal axes Rx and Ry of 
the corresponding eUipsoid by i?45 = [{Rl + i?^)/2]^/^. 
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FIG. 2: Comparison of the aspect ratio Ri^/L of a par- 
tially crystallized ^Be"*" plasma and the cold fiuid prediction 
for zero-temperature ellipsoidal plasmas, as a function of the 
ratio u^/uy between the axial and the smallest transverse 
trap frequency. Full squares: large crystals from FigQ Open 
(gray) circles: small (medium-size) crystals from Fig|3 Dot- 
ted line: J?45 / L for a low-density plasma in the gas phase with 
negligible particle interactions, leading to Rx/L = ujz/i^x, 
Ry/L — ujz/iOy. Inset: Cross-section of linear trap and of 
ellipsoidal plasma. Observation direction as in FigQ 

For the last crystal in the sequence, the boundary 
shows clear deviations from an ellipse; we attribute this 
to the presence of sympathetically cooled impurities 
of higher mass than ^Be"*", located at larger radii as 
compared to ^Be"*". The asymmetry along the trap axis 
is caused by cooling light pressure, which is not felt 
by the sympathetically cooled ions. In this case, the 
elliptical fit has been chosen to match the fragments of 
the outermost shell, which still exist at the left and right 
ends of the crystal. 

FiglH shows a comparison between the measured aspect 
ratio R45/L and the theoretical result from Eq.®. 
The agreement between experiment and theory is good, 
considering that the ^Be"*" ion crystals did exhibit two 
phases and were not pure. In addition, a systematic 
deviation between theory and experiment is expected 
at the largest applied voltages Vdc because then the 
smallest ellipsoid dimension, R^, becomes comparable 
to the shell spacing. In this limit, the cold fluid model 
continuum description is inaccurate. 
The cold fluid model also determines the relative change 
of the crystal length 2L when Vdc is changed. For the 
actual trap settings 2L is expected to decrease to 83% 
from its initial value, when changing Vdc from to 
4.2 V. The length of the crystal in the last image (e) in 
Figgis 84% of the initial one (a), in good agreement 
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with the expected compression. Therefore, while at 
large Vdc the observed transverse crystal shape starts 
to deviate from the predictions of the cold fluid model, 
the axial shape is still in good agreement, since the axial 
dimension of the crystal remains large compared to the 
shell spacing. 

While for large Coulomb crystals a spatially averaged 
description is a good first approximation, in small 
crystals one may expect effects related to the particle 
structure to show up clearly. As an example, FigOl (a)- 
(d) shows a crystal containing about 20 ^Be"*" ions and 
several sympathetically cooled low-mass impurities. In 
an effective trap potential with cylindrical symmetry, 
the crystal exhibits a sing le ^Be+ shell, FigEl(a). This 
shell appears smeared out in the CCD image (exposure 
time: 2 s), possibly because of a rotation- like diffusion 
of the ions around the trap axis. In the ellipsoidal crys- 
tals, Fig|3Ib)-(d), this diffusion is suppressed because 
it would require overcoming an energy barrier, and 
therefore the image shows individual ions. The most 
apparent property arising when the static voltage Vdc 
is increased is the completely dark region containing 
the sympathetically cooled particles. Furthermore, 
some of the fluorescence spots representing ^Be^ ions 
appear smeared out and show a reduced intensity. 
We attribute this to the micromotion occuring at the 
locations of these ions. As a test, we shifted the crystal 
by means of an additional static voltage applied to 
one trap electrode. This caused a rearrangement of 
the '^Bc^ relative to the dark core, but those ions that 
appeared well defined were always in the same region. 
Although the small crystal (a)-(d) in Fig|3| does not 
exhibit closed boundaries, we could still fit ellipses to 
them. In order to compare such small crystals with 
the cold fluid model, half of the typical shell spacing of 
29 yCtm was added to -R45 and L before calculating the 
ratio Ri^/L. The results are indicated by open circles in 
Fig|21 Even for this case there is a good agreement with 
the cold fluid model. Finally, in FigO (e)-(h) we also 
present a medium-sized crystal containing ~ 500 ^Be'*' 
ions and, again, additional low-mass impurities. The 
corresponding aspect ratio data as well as an additional 
data point for Vdc = 2.2 V, reported in Fig|21 also 
show good agreement with theory. However, for large 
and medium-sized crystals a deviation between theory 
and experiment becomes obvious as the values for the 
static voltage Vdc increase. Our molecular dynamics 
simulations show, that the observed deviations can be 
explained by small (< 15%) admixtures of sympatheti- 
cally cooled molecular impurities, in particular, /H^ 
ions close to the trap axis originating from residual 
gas contaminants, and BeH"*' ions located in the outer 
regions of the crystals formed by chemical reactions. 
The pronounced asymmetric ion distribution in the 
small, Figl^l (a)-(d), as well as in the medium-sized 
crystal, Figini(e)-(h), is not due to the broken cylindrical 



FIG. 3: Left: small crystal containing ~ 20 Be ions and 
a smaller number of sympathetically cooled impurity ions at 
different values of the static voltage VdC- V (a), 2.8 V (b), 
3.6 V (c), and 4.2 V (d). Right: medium-size crystal contain- 
ing ~ 500 ^Be+ ions. Vdc is set to OV (a), 1.4 V (b), 2.8 V 
(c), and 4.0 V (d). The asymmetric ion distribution in (b)-(d) 
and (f)-(h) is due to stray electric fields. 

symmetry of the trap potential, as this does not produce 
any visible asymmetry in the CCD images. Instead, 
we attribute the asymmetric ion distribution to stray 
potentials. While it was not possible to compensate 
for these imperfections by additional static voltages, it 
was always possible to reverse the asymmetry by means 
of these voltages. Our MD simulations confirm this 
interpretation. 

A direct estimate for the translational temperature of 
the Be+ is obtained from the spectral line shape of 
its fluorescence as the cooling laser is tuned towards 
resonance and the ion ensemble crystallizes. Since 
the temperature of the particles changes during the 
frequency scan, we fit a Voigt profile to each point of the 
recorded fluorescence curve to determine an upper limit 
for the Be"*" temperature. For small crystals (<1000 
particles), we find an upper limit for the temperature at 
the end of the scan of 42 mK. However, the accuracy of 
this method is limited due to the experimental resolu- 
tion. An indirect upper limit is obtained by comparing 
the size of the ion spots with MD simulations; here we 
find a tighter limit of <10 mK for the Be"'' temperature. 
We deduce, assuming thermal equilibrium, that the 
temperature of the sympathetically cooled impurity ions 
in FigElis <10 mK. 

In summary, we have studied the static behavior of 
Coulomb crystals in a fully anisotropic effective trap 
potential. We have found a good agreement with the 
simple cold fluid plasma model for small anisotropy. 
For larger anisotropy, deviations could be explained 
by the presence of additional, sympathetically cooled. 
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ion species. From an experimental point of view, the 
ability to reversibly deform a crystal permits to separate 
lower-mass sympathetically cooled ions from the laser- 
cooled ions. This allows to obtain a clearer picture of 
the impurity ion ensemble, without any background or 
foreground fluorescence from laser-cooled ions. This will 
also permit to manipulate the sympathetically cooled 
ions in a more direct way. The ability to generate a 
variety of ellipsoidal crystals opens up several directions 
for further study, e.g. oscillation modes of such crystals, 
and in particular the modes of two-species crystals. 
These modes could be of importance for identification 
of the non-fluorescent species. On the theoretical side, 
it is of interest to perform detailed studies of structures 
using MD simulations, which are able to take fully 
into account the particle nature of the cold plasmas. 
As an example, our MD simulations have shown that 
in strongly squeezed ellipsoids closed ion rings can 
occur. These represent a novel form of low-dimensional 
artificial structure whose detailed investigation should 
be of significant interest. 

We thank H. Wenz for the MD simulations. This work 
was supported by the Deutsche Forschungsgemeinschaft 
and the EU Network "Ultracold Molecules". 
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